Diamagnetism of nodal fermions 
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Free nodal fermionic excitations are simple but interesting examples of fermionic quantum criti- 
cality in which the dynamic critical exponent z — 1, and the quasiparticles are well defined. They 
arise in a number of physical contexts. We derive the scaling form of the diamagnetic susceptibility, 
X, at finite temperatures and for finite chemical potential. From measurements in graphene, or in 
Bii_xSbx {x — 0.04), one may be able to infer the striking Landau diamagnetic susceptibility of 
the system at the quantum critical point. Although the quasiparticles in the mean field description 
of the proposed d-density wave (DDW) condensate in high temperature superconductors is another 
example of nodal quasiparticles, the crossover from the high temperature behavior to the quantum 
critical behavior takes place at a far lower temperature due to the reduction of the velocity scale 
from the fermi velocity vf in graphene to y/VFVoow, where wddw is the velocity in the direction 
orthogonal to the nodal direction at the Fermi point of the spectra of the DDW condensate. 
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I. INTRODUCTION 

In a class of quantum critical points (QCP), Lorentz 
invariance appears as an emergent symmetry, but in gen- 
eral the quasiparticle residue, as inferred from the one- 
particle Green's function may vanish. In rare cases, when 
the quasiparticle residue is finite, depending on the statis- 
tics of the excitations, the Lorentz invariant QCP is de- 
scribed by either a relativistic massless bosonic free field 
theory (massless Klein Gordon action) or a relativistic 
massless fermionic free field theory (massless Dirac ac- 
tion). Only in (l+l)-dimcnsion are both descriptions 
identical due to transmutation of statistics. Though the 
theory has a relativistic form, the speed of excitations is 
usually about two orders of magnitude smaller than the 
physical speed of light. Due to fluctuations on all length 
scales in a critical system, many physical quantities ex- 
hibit power laws and obey scaling in the vicinity of the 
QCP. Even in the simplest of such systems, there are 
surprises buried in their diamagnetic response because a 
magnetic field is never a small perturbation: any per- 
turbation that changes the spectra from continuous to 
discrete can not be considered small. Here, we hope to 
elaborate on this topic and present estimates that may 
be tested in experiments. 

For a class of tight binding models in the half filled 
limit, for example graphite or graphene, the energy van- 
ishes at distinct points of the Brillouin zone known as 
the nodal pointsjii^ and in the long wavelength and low 
frequency limit the dynamics are well described by Dirac 
fermions obtained by linearizing the spectrum around the 
nodes. The nodal spectra can also arise from a conden- 
sate. An example is nodal fermionic quasiparticles of 
a particle-hole condensate in ^ ~ 2 angular momentum 
channel, as in a singlet d-density wave (DDW), staggered 
flux phase, or an orbital antiferromagnetf2ili^ 

The electromagnetic charge is a conserved quantity for 
a tight binding model of an electron. This is also true 
if the order parameter is a particle-hole condensate, as 
in a DDW. In these cases, the electromagnetic field can 
be incorporated via the minimal gauge coupling. We 



shall restict ourselves to such systems and not consider 
nodal Bogoliubov quasiparticles of a d-wave supercon- 
ductor. The contrasting response of d-wave superconduc- 
tor (DSC) and DDW is evident.^ The quasiparticles in 
a superconductor do not minimally couple to the vector 
potential A, but to the supercurrent ^ {Vip — 2eA/hc), 
where Lp is the phase of the superconducting order pa- 
rameter, e is the electronic charge, and c is the velocity 
of light. 

The effect of the chemical potential, /i, is extremely im- 
portant, as it can introduce electron or hole pockets and 
render the linearized free Dirac theory invalid. However, 
for small /i one can still use the linearized continuum the- 
ory; /i = describes the vaccum of the relativistic mass- 
less theory and hence is critical. But, for a finite /z, one 
is dealing with a finite density of excitations. Thus, one 
is perturbed away from the criticality, and this should 
provide a cutoff. 

For the diamagnetic response at /i = and zero tem- 
perature (T = 0), one can use a simple quantum critical 
scaling analysis to find the power laws satisfied by the 
magnetization and the susceptibility.^ In this paragraph 
we shall set e — fi, = c = 1. From gauge invariance, the 
vector potential A has the same scaling dimension as the 
momentum which is where L is a length. Therefore 
the magnetic field H has the scaling dimension or 
there is a length scale L ~ H~^^^. One can immediately 
see that this length, which acts as a cutoff at the quantum 
critical point of the free Dirac fermions, is proportional 
to the Landau length. Since the hyperscaling should be 
valid for d — 2 and the dynamic critical exponent z = 1, 
the singular part of the ground state energy density, fio, 
multiplied by the correlation volume, L^'^^^\ should be 
a universal number;^ that is, Qo ^ H^^^. Therefore, the 
magnetization behaves as M ~ — iJa and the diamag- 
netic susceptibility behaves as x —H~2. In the H ^ 
limit, X diverges, which will be cut off by a number of 
physical effects not contained in this argument, and the 
stability of the state may not be in question. 

The diamagnetic sign cannot be obtained from the 
scaling argument. The energy levels in a magnetic field 
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are bunched (discrete spectra) although the mean den- 
sity of states is unchanged. The number of quasiparticles 
that can be accommodated below any given energy de- 
pends on whether or not this energy coincides with an 
eigenvalue of the Landau spectrum or falls in between 
two eigenvalues. For the nonrelativistic case it is easy to 
see that on average the energy is increased, because near 
E = we always start with an empty interval. For the 
relativistic continuum theory of nodal fermions where we 
have to impose an ultraviolet cutoff, this is subtle and re- 
quires a proper regularization. Using the work of many 
authors involving ^-function regularizatioufSwe can show 
that the answers are indeed cutoff independent and the 
energy is increased. This was also checked by considering 
a lattice version and Peierls substitution to incorporate 
the magnetic field. 

We cannot apply the same scaling argument to free 
Dirac fermions in (3 -I- l)-dimcnsions because hyperscal- 
ing is violated. This case is best described by a mean 
field theory with logarithmic corrections. It is known 
from explicit calculations that the singular part of the 
ground state energy density fip ~ H^logH.^^ A naive 
application of the above scaling argument gives only a 
regular contribution, flo ~ H^, which is not surprising. 
Thus, we feel confident that the quantum critical scaling 
analyses are indeed meaningful. 

Consider d — 2; some aspects of the finite temperature 
and finite chemical potential results can be understood 
from the notion of quantum criticality. From finite size 
scaling, the correlation length, ^(T), is proportional to 
the thermal wave length, 

Xt = (1) 



that is. 



hvF 



(2) 



where Aq is a universal number of the order of unity and 
is the Fermi velocity. Tuned to ^ = 0, the quantum 
criticality will persist until £,{T) is the order of the lat- 
tice spacing a. Since vp is large, one would naively expect 
the singular diamagnetic susceptibility X cx -iJ-V2 to 
persist over a wide temperature range. In fact, x is gov- 
erned by a balance between two length scales: the Lan- 
dau length, Ib = (Sc/2eS)i/2, and ^(T). If ^(T) > 
then X follows the power law indicating the quantum crit- 
ical behavior, and in the opposite limit we obtain linear 
response, x ^ —1/T. At T = 0, non-zero fi tunes the 
system away from criticality. For small fi one can still 
use the linearized spectrum, and this introduces another 
length scale, which is essentially the interparticle spacing 
A ~ hvp/lJ-- For X> Ib, X follows a power law and in the 
opposite limit x ^ "I/a*- 

Quantum criticality of relativistic Fermions is exper- 
imentally relevant for graphene for which a linear spec- 
trum has been established experimentally."'^^ These quasi- 
particles are charged fermions and show anomalous inte- 
ger quantum Hall effect, as well as Shubnikov-de Haas 



oscillations J^ii^ii^ It is then natural to expect that, as 
T ^ 0, graphene should have the signature of a dia- 
magnetic "instability" consistent with quantum critical- 
ity described above. Similarly, the diamagnetic suscepti- 
bility of Bii_xSbx {x = 0.04) for which the linear disper- 
sion of the fermionic excitations is known to be presenti^ 
remains unexplored. This should be approximately de- 
scribable in terms of a (2 + l)-dimensional Dirac theory 
with weak interlayer couplingiii As mentioned above, it 
has been suggested that the pseudogap phase of the high 
Tc superconductors can be described by DDW, whose 
quasiparticle excitations for fi = are Dirac fermions, as 
was recognized a long agoi^ Our work is an extension of 
these early analyses of diamagnetism of nodal fermions to 
finite temperatures and finite chemical potential, which 
leads to interesting results. 

In a set of magnetization measurements, Ong and his 
collaborators-'^* have uncovered unusual diamagnetism in 
the pseudogap state of the high temperature supercon- 
ductor Bi2Sr2CaCu208+x (BSCCO). In the pseudogap 
regime, above the superconducting transition temper- 
ature, the diamagnetic susceptibility diverges as x ~ 
_jj{i-S)/s^ 0, where the effective exponent (S(T) 

is greater than unity over a very broad range of tem- 
perature. Such a divergent susceptibility above a phase 
transition calls for new ideas, because the response in 
general should be linear. Only at a critical point, where 
there are fluctuations on all scales, is it possible to ob- 
tain such a nonlinearity. In particular, it is known that 
for two-dimensional Kosterlitz-Thouless theory S — 15 a,t 
criticality;^ T = Tkt, but the response is linear for any 
temperature T > Tkt- To the extent that the critical re- 
gion is sufficiently wide, it is of course possible to obtain 
a large value of susceptibility, but not a divergent suscep- 
tibility, as seen in measurements where fields as small as 
5 Gauss were used. Taken at its face value, experiments 
indicate a critical phase extending over a wide region of 
the pseudogap state. 

Long ago it was suggested that a weakly coupled stack 
of XF-systems could exhibit a floating phase in which 
the three-dimensional behavior at low temperatures con- 
verts to a floating power-law phase (a stack of decoupled 
layers) at intermediate temperatures and flnally to the 
disordered phase at high temperaturesi^ It is now rig- 
orously known— that if the coupling between the layers 
is Josephson-like (a likely scenario), a floating phase is 
ruled out even for arbitrarily long-range couplings. Very 
special, finely tuned, interlayer couplings are necessary to 
produce a floating phase, which appears to be unlikely. 

Although we flnd that a sizable diamagnetism sets in 
with the DDW gap over and above the conduction elec- 
tron diamagnetism, our results cannot explain the data of 
Ong and his coworkers: (a) there is no finite temperature 
critical phase; (b) the relevant scales are vastly different. 
As mentioned above, Kosterlitz-Thouless theory cannot 
account for a critical phase above Tc, though the order of 
magnitude is reasonably close. We hope that our calcu- 
lated crossover behavior of the diamagnetic response will 
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be observable, at least in graphene or in Bii_xSbx- 

The paper is organized as follow: in Sec. II we will 
describe the effective model for nodal fermions in two di- 
mensions and outline the formalism for computing the 
grand thermodynamic potential. In Sec. Ill wc will de- 
scribe our results for two dimensions (2Z?). We first de- 
scribe the results for the case /i = and then proceed 
to the discussion of 7^ 0. In Sec. IV we consider weak 
interlayer coupling in the context of a three dimensional 
(3-D) system. In Sec. V we consider numerical estimates 
of the effects that are experimentally relevant and in Sec. 
VI we conclude. There are two appendices that contain 
certain mathematical details. 



where Dj = dj — i^Aj is the covariant derivative. Landau 
levels can be easily found by squaring the hamiltonian to 
be 

En = zb^VH = ±VaBn (7) 
Ib 

where Ib — (Sc/2ei3)^/^ is the magnetic length. We 
have introduced a = 2hevp/c for notational clarity. The 
same formalism can be applied to the nodal spectra of 
Bii_xSbx {x = 0.04). 

B. d-density wave 



II. NODAL FERMIONS: TWO-DIMENSIONAL 
SYSTEMS 

A. Graphene 

When linearized about the two inequivalent vertices 
of the Brillouin zone, the tight binding Hamiltonian, H, 
defined on a honeycomb lattice of a sheet of graphene 
involving only nearest neighbor hopping, with matrix el- 
ement t, becomes in the continuum limit (lattice spacing 
a ~* such that at is finite) 



H 



-hvf. 



(Pk 



s2, 



(3) 



where ipi and ■02 are two species of two-component Dirac 
fermions corresponding to two inequivalent nodes, and 
vp = V3at/2h is the Fermi velocity; the spin index s is 
summed over. The sum over two inequivalent nodes can 
be written in a compact and Lorentz invariant form as 



H = —ihvi 



2 



(4) 



where ip = ip'^'^'^ and ip — 



is now a four component 



spinor, ignoring the irrelevant spin indices. We are using 
a reducible representation of 7-matrices formed from the 
standard Pauli-matrices cr's: 
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-i(Ji 



,7' = 



1(72 









(5) 

The Landau level problem in the tight binding formu- 
lation is a Hofstadter problemi^ But, for weak enough 
magnetic fields we can analyze the continuum model by 
incorporating the magnetic field by minimal coupling pre- 
scription. So, the hamiltonian of interest takes the form 



H 



-ihv, 



2 

•E 



d^xip^-' Djip 



(6) 



The nodal spectra of the DDW is also a well studied 
problem.— 1^ The low-energy quasiparticle Hamiltonian for 
the DDW state is 



H 



DDW 



d^k 

(2^ 



[(e(fc)-Ai)c^^(fc)cs(fc) + 

tWiky\k)csik + Q)l (8) 



where e(fc) is the single-particle energy, commonly chosen 
to be 

e(fc) = — 2i (cos k^a + cos kya) + At' cos k^acos kyU, (9) 

and Q = (7r/a,7r/a). The nearest neighbor hopping ma- 
trix element is t and the next nearest neighbor matrix 
element is t' . The spin-singlet DDW order parameter 
takes the form: 

{c'^(k + Q,t)cs'{Kt)) =iW{k) SI,, (10) 

where the gap function is given by 

Wo(T) 

W{k) = — - — (coskxtt — coskyo). (11) 

As the order parameter breaks traslational invariance 
by a lattice spacing a, it is convenient to halve the 
Brillouin zone and form a two-component Dirac spinor. 
Then, in the reduced Brillouin zone, the mean field 
Hamiltonian is 



H = 



d^k 



i(6(fc) + 6(fc + Q))-/x 



-(e{k) - e{k + Q))a3 + W{k)ai 



where 



X2s 



Cs{k) 
ics{k + Q) 



Xs{k), (12) 



(13) 



The spin index s can again be dropped, as this will not 
enter in our calculation except for an overall multiplica- 
tive factor. 

The quasiparticle energies are 



E±{k) 



1 



{e{k) + e{k + Q)) 



4 



e{k + Q)y +AW^ik). (14) 



At half-filling, ii — 0, there are 4 gapless nodal points 
at (±^,±^), the Dirac points. A non-zero value of 
fi will open up fermi pockets. The low-energy physics 
will be dominated by these gapless fermionic excitations. 
We choose a single pair of nodal points, (5^,5^) and 
(—52;,—^) and include the other pair of nodes into our 
final result. We take the x-axis to be perpendicular to 
the free-electron Fermi surface and the y-axis parallel to 
it at one antipodal pair of nodes; similarly, the x-axis is 
parallel to the free-electron Fermi surface and the y-axis 
is perpendicular to it at the other pair. Linearizing the 
spectrum about the nodes, the dispersion relation is 



E{k) = ±h^vlkl + v 



2 1.2 

DDW'^'i/ 



(15) 



where vp — 2^/2ta/h and wddw = Wa{T = 0)a/V2h. It 
is important to note that the parameter t' does not enter 
at linear order. It is now obvious that the formalism 
is identical to that described in the previous subsection 
provided we replace vp by (w_fWddw)^^^ and rescale k^ —> 
kx \/vuDW /vp and ky — + ky\Jvp j wddw to account for 
the DDW gap anisotropy. 



C. Grand canonical potential 

Consider the grand canonical thermodynamic poten- 
tial per unit area of a two-dimensional (2D) system: 



n(T, [i) = -fcsT / deD{e) log 2 cosh 



IkBT 



(16) 



Here, Die) is the density of states (DOS), which in the 
presence of an applied perpendicular magnetic field B, 
takes the following form: 



D{e) = CB 



(17) 



where C — Nfe/hc is an universal constant, such that, 
CB represents the Landau level (LL) degeneracy factor, 
i.e. the magnetic flux per unit area due to the applied 
field measured in the unit of flux quantum. Nf is the 
number of electron flavors - Nf — 4 for both graphene 
and DDW. Note that in Eq. (|17p we have assumed a pure 
system. The presence of disorder broadens the sharp 5- 
functions in D{e), however, we restrict our discussions to 
a clean system in this paper for simplicity. 

Substituting D{e) in Eq. (fTB|) . we can write n(p,T) = 
f2o(/x) + flT{lJ-), where ^oil^) is the temperature indepen- 
dent part (hence contributes even at T = 0) given by 



^o(a^) 
CB 



n— Tic + l 

-f^inc + ^) - V^Bi/2^(-i, 1 -I- nj- (18) 



Here we assumed fj, > (electron doping), and thus the 



positive LL's are filled only up to 



lnt[^J.^/aB] at 



T = while all the negative LL's are filled (Int[-] stands 
for the 'integer part'). Here C{s,q) — J2'kLo(^ + 
the standard Hurwitz ^-function. It is straightforward to 
/i < 0. The T-dependent contribution is 



»T(/i) 

CBkeT 



log(l + e ''b'^ 



log(l + e "bT 



oc 

log(l + e" 

n— ric + l 



e ''bt ) 



(19) 



Note that at finite T the thermal energy can excite elec- 
trons across /it to arbitrarily high (positive) LLs, and thus 
the n-sum must include the whole of Dirac cone, as shown 
explicitly in Eq. 



III. RESULTS: TWO DIMENSIONS 
A. Undoped system, /i = 

Conside the half-filled system: /i = 0, hence ric — 
At any temperature, the length scale of the critical fiuc- 
tuations is the correlation length £,{T). Thus, in order 
to observe the T — Q critical behavior the largest length 
scale for the system must be this length. In the pres- 
ence of a magnetic field B, the response of the system 
will show critical behavior only when ^(T) > Ib- At 
T = 0, this condition is trivially satisfied, because the 
length scale of the critical fiuctuations is infinite, and we 
obtain 



1 

CaB^^^i^ 



71=1 

CySB3/2<;(3/2) 



-A/"o7V/g2DA'lv^S3/2. 



(20) 

(21) 
(22) 
(23) 



Here ({s) — J2 



k=l ' 



is the Riemann (^-function, Aq = 
3C(3/2)/87r « 0.312, = eh/2mc, the Bohr magneton, 

•'Vp/HB 



with m the free electron mass. The scale Bq = 
is a material dependent constant and has the dimension 
of a magnetic field. The transition from the second ex- 
pression to the third is an example of standard (^-function 
regularization of a divergent sum over n. The proof fol- 
lows from the remarkable result due to Riemannf2i that 



2^-"r(s)C(s)cos ( isTT 



(24) 



The logic is that the "divergent" sum is physically cut 
off at some value of n and is not truly divergent, but 
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a gauge invariant regularization is necessary. This is ac- 
complished by the analytic continuation given by the Rie- 
mann reflection in principle. Other regularizations are 
given in Refs. andfiol. 

For reasons of physical transparency we shall often 
express our formulas in terms of an equivalent non- 
relativistic free electron gas, while keeping in mind that 
the real parameters that enter our calculations, such as 
vf, Nf, etc. bear no real relation to this free electron 
system with a circular Fermi surface and two flavors of 
spin. Thus, we have written 



(25) 



which is the standard, energy independent DOS of a two- 
dimensional {2D) non-relativistic Fermi gas. Similarly, 
we can express 



/52D = g2Dmvp/2 



(26) 



where p2B is the 2D arial density. Here we have used the 
transcription vp = tikp/m, where i kp is the Fermi wave 
vector of the equivalent non-relativistic Fermi gas. 

The corresponding T 7^ contribution takes the fol- 
lowing form: 



log2-|-^log(l + e 



(27) 



It becomes clear from Eq. (|27|) that fix is a function of 
the ratio of the two fundamental length scales Xt/Ib, and 
thus it must have a scaling form. 

We calculate the magnetization M and the suscepti- 
bility X from 



M 
X 



dM/dH, 



(28) 
(29) 



where H is the magnetic field strength. These also have 
scaling forms. If we introduce 



b = Xt/Ib, (30) 



we obtain: 



(31) 



X = Xo + XT 

3CVSC(3/2) 

Xo =- 



XT = Xo^^^{S+^idS/db)} = xofib), (33) 



where S is given by 



n=l 



_j_ g6\/rr 



(34) 



The function f(b) defined in Eq. (|33p is a universal func- 
tion of its dimensionless argument, which can be written 



as a series expansion in b (See Appendix A): 



fib) 



An 



am 



b_ 8 

18 ~ 



t (4g + i) 



X {q+l)rj{-Aq-l)a-2q-l) 



(35) 



where, 77(5) = [1 — 2 ]C(s) is standard Dirichlet 77- 
function. 

In the limit S> Is, or equivalently 6 3> 1, it is 
the quantum criticality that dictates the response of the 
system, and Eq. p5p is not particularly useful. Instead, 
we can obtain the analytic expression for f{b) in this 
regime by replacing S = 2 J2n ^^Pi^b\^) in Eq. ((34|) 
to get (See Appendix B for details) 



fib) = Fy^ib) - b^F,^,{b) + -Fy,{b), 



(36) 



where we have defined a (convergent) 6-dependent inte- 
gral Fp{b) as, 



Fpib) 



4V¥ 
C(3/2) Jo 



dx- 



- 1 



(37) 



In fact, it is possible to obtain an explicit 6-dependence of 
X by estimating the saddle-point approximation of -Fp's, 
which results in 



x(At > Ib) ~ Xo 



1 



N. 



3/2 



N, 



5/2 



b \ (e''V6 _ 1) 



(ebViO-i) 12(e''Vi4_i) 
where, Np is a pure constant given by 



" C(3/2) 




^-)9[l-log(4p)] 



(38) 



(39) 



The message from Eq. (|38|) is transparent: for & 3> 1 the 
first term, xo: dominates, causing the B~^^^ behavior in 
the susceptibility, while the rest of the terms in xt van- 
ish exponentially. As 6 is decreased, xt grows, modifying 
the non-linearity of x in B. This behavior continues until 
6-1, that is, until ^(T) ^ Xt ^ Ib- Finally, for 6 < 1 
the critical fluctuations fail to describe the magnetic re- 
sponse, and the susceptibility follows linear response. For 
6^1, can we keep only the first two terms in Eq. (|35p : 
the next term is ~ b^ and hence negligibly small. The 
first term exactly cancels xo and we have. 



X = 



Ca 
2Ak^T 



= -N 



P2T>l4} 

ikBT ' 



(40) 



It is the expected diamagnetic, i?-independent behavior 
in the high temperature limit if we absorb Nf in the 
definition of the areal densityi^ 

We plot log/ (6) as a function of log 6 in Fig. 1, using 
three following methods: (a) by numerically evaluating 



6 



XI -5 



6jO 
O 



-10 



o large b approximation 
• small b approximation 




Ib > Af 



4.10.26.27 



FIG. 1: log[— /(&)] as a function of log 6. The numerical 
evaluation of Eq (I34p is given by solid line, and the analytical 
expression for the large and small 6 limit are given by empty 
and solid circles respectively. The inset shows the blown- 
up crossover region (in linear scale). The two asymptotic 
limits reproduce f{b) surprisingly well over almost the entire 
parameter regime. 



Eq. ([34]) with a desired (high) accuracy for a wide range 
of b; (b) from the large b asymptotic expression as in 
Eq. ([311), and (c) evaluating Eq. dSS]) in the hmit 6 < 1, 
which amounts to keeping only terms up to linear or- 
der in 6. We find that the two asymptotic expressions 
encompass almost the entire parameter space surpris- 
ingly well. The smoothness of f{b) implies that while 
X ^ B~^/^ for 5^1, its behavior smoothly crosses over 
to i3-independent diamagnetic behavior for 6 < 1. Note, 
however, that f{b) itself is finite at all b, and thus the 
Landau diamagnetism prevails. 



B. Effect of finite 

When the doping is small, fi is small as well and cor- 
responds to an effective quasiparticle description. First, 
consider T = 0; as long as /x is small enough for the lin- 
earization of the spectrum to be valid, the Fermi surface 
changes from a point in momentum space for /i = to 
a circle, and generates a length scale of A = hvp/fJ-, the 
inter electronic spacing. In the limit X> Ib we get from 
Eq. HID 



Xo(Ai) 



1 



(B + ^2/a)i/2 



(41) 



in the leading order. It is now obvious that for X ^ Ib 
we get X B^^/^. This divergence of x is cut off for 



and we get 

Ca 



(42) 



This is of course expected because the chemical potential 
tunes the system away from the quantum criticality. For 
finite B, in the (non-critical) regime of X < Ib we expect 
de-Haas Van Alphen (dHVA) oscillation in the magneti- 
zatioE^i^iSl due to the cutoff introduced by /x. 

For T ^ 0, the additional T-dependent part in Eq. (fT9|) 
becomes important; see Eq. (jASp . Because we now have 
three different length scales: Xt, Ib and A, the expression 
for (and x) will depend on their relative magnitudes. 
The most important regime from the perspective of crit- 
icality, <C A <C At, is particularly simple. In this 
case we can use similar approximations as in Eq. 
yielding: 



l + i<{ cosh(AT/A) 



Xih < A < At) = xo 
X (^Fy2ib)-b'F,/^{b) + ^Fj^,ib)^ 



(43) 



Thus, the susceptibility has a scaling form in terms of 
two independent dimensionless variables: At/ A and b = 
Xt/Ib- The expression for x i^i Eq. ()43|) is valid even if 
<C At < A, but the latter condition invalidates the 
applicability of the linearized theory due to large /x. 

In the opposite limit of linear response, simple expres- 
sions for the susceptibility can be derived, and we get 



(44) 



which reducuces to Eq. ((40|) when /i = 0, as expected. 



IV. THREE DIMENSIONS: EFFECT OF WEAK 
INTERLAYER COUPLING 

Materials where this 2D nodal fermion theory is appli- 
cable are layered (quasi-2D) systems, an exception being 
graphene, which is indeed atomically thin. If wc include 
weak interlayer coupling in a tight binding hamiltonian, 
the energy spectrum acquires an additional quadratic dis- 
persion given by 



E{k) = tx_klf ±hvFJkl 



1,2 



(45) 



where t±_ is the interlayer hopping matrix element and 
£ is the interlayer spacing. Introduction of this new en- 
ergy scale will cut the divergence off x(r ^ 0) when 
the magnetic energy scale becomes smaller than t±_ . The 
corresponding Landau energy spectrum is 



(46) 
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For T = 0, /.t = limit, we get 



CB 



dk^ V k_Lfc^^2 _ ^aBn (47) 



where ric = Int 
we get 



{t^kie/^aBY 



Performing the n sum 



CB 



tt/I 



t^kiec{Q,fi, + i) 



iBC(-l/2,n, + l) 



(48) 



If [t^kli"^ l^/aBY < 1 for any value of k^ within the 
cutoff, he = 0, and the kz integrals can be done trivially. 
Thus, 



3D 



CBtiH^ CySB3/2C(-l/2) 



I 



(49) 



The susceptibility now is just the previous zero tempera- 
ture result divided by This implies hvp/lBtA. < ti"^, 
which leads to a lower cutoff in the magnetic field, 
given by B^ = TT'^t'^c/{2ehvp). For a given t±, and 
B > Be, X ^ B^^l^ . When is vanishingly small. 
Be is also vanishingly small and can be ignored. When 
\_k?z^ I \foiBY > 1 for any value of kz, the result is 
more complicated and will be representative of a truly 
3D system. 

However, in 3D electrodynamics one has to distinguish 
between B and H, which leads to another cutoff. Follow- 
ing Ref. 0, we provide the appropriate formulas for x 
/i = 0. In 3D electrodynamics the magnetic induction B 
and H must be distinguished: 



B = H + A7:M3d{B). 



(50) 



For X we must find _B as a function of H. Since, M in 
general is a function of B and T, _B is a function of H 
and T. From Eq. ((231) we get for Xt :> Ib 



M3D{T = 0) = -2J\foNfg2BfJ-lVB^/i (51) 

Now using Eq. ((5T|) we obtain the relation between B and 
H: 



B{H,T = {))^[{H + H,Y/^-Hl'^f (52) 



where, 



(53) 



and has the dimension of H. Plugging Eq. ([5^ into 
Eq. jnil) we get. 



-1/2 



(54) 



The same analysis in the linear response regime, At ^ 
Ib, yields 



and 



X3d(At < Ib) 



H 



1 + To/T 



1 1 



(55) 



(56) 



where Tq = AiiNfp2Dii^g/3kBi- 



V. EXPERIMENTAL RELEVANCE 

We have established in Sec. IIIII that the diamag- 
netic susceptibility undergoes a crossover as a func- 
tion of T, from its zero temperature power-law behav- 
ior {x ^ —B~^/'^) to high temperature linear behavior 
(x ~ ~l/2^)- It is interesting to ask if this crossover is 
observable. Consider graphene; we take the experimental 
value of vp — 10^ cm/s and use Eq. ([^0]) for the high-T 
regime. We obtain 



X2D 



9.88 X 10 



-10 



/ 2 

emu/ cm 



(57) 



where the temperature T must be expressed in Kelvin. 
Note that X2D = X of Eq. 

In order to compare with experiments on layered (quasi 
2D) materials we calculate the susceptibility by dividing 
Eq. [301 by L If we now take I = 3.35 A (the value for 
graphite), we obtain 



X2D 

I 



2.95 X 10" 
T 



The corresponding susceptibility per unit mass is 
X2D 0.0134 



T 



-cmu/gm. 



(58) 



(59) 



using the mass density of graphite^ — 2.22 g/cm^, 
which agrees very well with the experimental results.— 

Upon lowering T, X2D is strongly enhanced and the 
power-law region can be accessed for At > Ib- This 
implies that in graphene for B > 5.6 x 10~^T^ (Gauss) we 
must use Eq. (|38p instead of Eq. ([^0]) for the estimation 
of X2D- In particular, in the T — > limit we obtain 



X2D 



3D 



0.012 



emu/gm. 



(60) 



We demonstrate in Fig. 2 the behavior of X2Ti/^P^d ^s 
a function of T by numerically evaluating Eq. ([3^ to 
illustrate the aforementioned crossover behavior. 

However, if we use Eq. ([M)) to take into account the 
demagnetization effect due to interlayer coupling in 3D 
graphite, we obtain 



X3D ~ - 



2.95 X 10-2 
(T + 0.37) 



(61) 
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FIG. 2: Evolution of x = X2d/^P3d (defined in text) 
as a function of T for various values of _B. It is calcu- 
lated for graphene, using vf = 10* cm/s, £ = 3.35A and 
pTd = 2.22gm/cm^ 



and this is a very small effect for high temperatures. For 
3D graphite, when At > ^s, we use Eq. ([53]) and Eq. ([54]) . 
For graphite iJ, =0.11 G and 



X3D = 



47r (1 + 9.m)i/2- 



(62) 



Therefore, in the limit of <C 0.11 Gauss graphite 
should become a perfect diamagnet (!), which, however, 
is a very small field. If the condition, Xt > Ib com- 
bined with the value of the scale H^, we find that when 
T < 1.5K, the demagnetization effect will be important. 

In Sec. Ill Bl we described the DDW phase of high Tc 
superconductors by the Hartree-Fock theory of the nodal 
fermions in the copper-oxide layers. We shall estimate 
the strength of the diamagnetic susceptibility from the 
DDW order, using the following experimental parameters 



for typical cuprates 



18,30 



3 X 10^ cm/s and £ = 12A, 



where wddw is estimated assuming a fully formed DDW 
gap Wo ~ 35 meV which leads to the anisotropy in the 
velocity i'.f/wddw ~ 28.6 if t = 250 meV. 

In the linear response (x ~ regime, we obtain 

from Eq. ((M)). 



2.6 X 10^5 
'(r-F3.3 X 10' 



2.6 X lO"'^ 



When Xt > Ib, we use Eq. 
= 2.7 X 10-5 G and 



1 



X3D 



(63) 



53]) and Eq. ^ to obtain 



(64) 



47r(l + 3.7x 104i/)i/2 



This indicates that the diamagnetic susceptibility of 
DDW from nodal fermions may be measurable. The dis- 
cussion above would imply that DDW would become a 
perfect diamagnet when H <C 2.7 x 10"^ G, which, how- 
ever, is such a small field that many other effects will 
intervene, and one would observe x ~ —H^^^^, but not 
perfect diamagnetism. 

In the experiment on BSCCOi^ a T-dependent power 
law is observed over a wide range of temperature in the 
small H limit. Moreover, at the smallest value of the 
magnetic field in Ref. [l^, = 5 G, we get from the 
DDW calculations X3D « — 1.9 x 10"'' in CGS units as 
T — > 0. This is orders of magnitude smaller than that 
found in the experiment. Therefore, the magnitude of 
the diamagnetic susceptibility of Ref. Ts' can not be ex- 
plained within a DDW framework alone. One must note 
however that as the temperature is lowered, the system 
will generically enter from the DDW phase to a coexist- 
ing DDW and d-wave superconducting phase for much of 
the parameter regimci^ Thus it is clear that the super- 
conducting diamagnetic effects of the Kosterlitz-Thouless 
theory cannot be ignored.^- But of course none of these 
considerations can explain the observed critical phase, 
which requires new ideas. 



VI. CONCLUSION 

We have shown that the notion of quantum criticality, 
although restricted to non-interacting nodal fermions as 
elementary excitations, offers interesting insights to dia- 
magnetism of semimetals. When the the chemical po- 
tential is zero, the system is inherently quantum critical, 
and we derived the scaling function for The scaling 
form suggests that the non-linear behavior of x as a func- 
tion of B, due to quantum criticality, can persist up to 
a large enough temperature, which may be accessible in 
measurements in graphene. We have also discussed how fi 
tunes the system away from the quantum critical region. 
The root of the large magnitude of the diamagnetic sus- 
ceptibility in graphene or graphite is of course the large 
Fermi velocity vp- 

There are a number of difficult but obvious questions 
regarding the roles of electron-electron interaction and 
disorder. These could be topics for future work. We 
have seen that our simple picture of the DDW does not 
explain the remarkable experiments in the high temper- 
ature superconductors. We do not know if the general- 
ization of the Hartree-Fock picture of the DDW to the 
six-vertex model where a power law high temperature 
phase was found^^ will be able to explain these experi- 
ments. It is certainly worth exploring. We stress, for the 
reasons stated above, that these experiments are not fully 
explained by Kosterlitz-Thouless theory, as is sometimes 
claimed. 

It is clear that the Euler-MacLaurin summation ap- 
proach to compute Landau diamagnetism for non- 
relativistic fermions fails because of the non-analyticity 
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due to massless Dirac fermions in semimetals. It is not 
known to us if there are any systems for which (3 + 1)- 
dimensional quantum critical behavior x ~ log-ff is ex- 
perimentally observable. The material Bii^xSbx is lamel- 
lar, as is bismuth telluride, and is better described as a 
two-dimensional system with weak interlayer coupling. 
Nonetheless, it would be interesting to study the dia- 
magnetism of this material as a function composition. 
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APPENDIX A: HIGH TEMPERATURE SERIES: 

We thank Stuart Brown and Subir Sachdev for discus- M 7^ 



The grand canonical thermodynamic potential is given by 



Vt{T,^i) = -CEksT 



K \ /J n— 1 ^ ^ ^ n— 1 ^ 



2kBT 



(Al) 



Each individual LL sum fails standard convergent tests. 
The technique to deal with such sums in the quantum 
critical regime is discussed in the text. The strategy in 
the other limit, where linear response holds, is to con- 
vert the LL sums to express them as series in powers of 
6 ~ l/T, so that meaningful conclusions could be drawn 
about the small h limit (equivalently, high T limit) by 
considering leading order terms systematically. The pro- 
cedure relies on C-function regularization, details of which 
could be found in literature, but the purpose of this ap- 
pendix is to provide a self contained description. Sep- 
arating out the zero temperature part f2o(M) ^md finite 
temperature part ilrifJ-) we obtain Eq. and Eq. (fT5)) 
respectively. We now wish to express JItCm) as a series 
expansion in powers of b. For this purpose we focus below 
to one term in Eq. say, the following one: 



Ii and I2 respectively), first the logarithms in powers of 
the exponentials and subsequently e t in a power 
series to write 



n=l 



00 00 , _.j-i£_,t 00 ,j. 



k — ' r! 

n— 1 k—1 r— 

OD ,„ 00 . 00 

^ (— e ''bt 
71" ^ 

r=0 ■ k=l 



^-Lii_,(-e-^)C(-r/2). 



(A3) 



r=0 



/ = ^ log(l + e ''"bt ) 

00 00 

= ^log(l + e"'^)- log(l 



Jic + l 



(A2) 



We will now expand both the summations (we call them 



In the third step above, we interchanged the order of 
summation, which in general leads to a correction, but 
in this particular case it is zero (for details see Ref. [90 . 
And in the final step we have used the standard definition 
of Polylogarithm Lis(z) = X]fc°=o ■^'^/^^ and the Riemann 
(■-function. Similar manipulations for I1 lead to the fol- 
lowing: 



=0 



00 

T 

k=l 



E-Lii-r(-e ^BT)(^(-r/2,l + n,) 



(A4) 



where we get Hurwitz's (^-function instead of Riemann (-function. Employing similar simplification to each term of 
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Eq. (fT9l) . we finally obtain the desired high temperature series expansion: 
nrifi) = -CBkeT log (l + e"^) - ^ -Lii_,( 
X {Lii_,(-e"^) + (-irLii_,(-e^)} 



)[i + (-inc(-^)+E7C(-^,i+"c) 

r=0 



We arrive at Eq. (|35p by letting = 0. Also, starting After carrying out the r-sum we get 
from Eq. (jASP we can derive Eq. (|^ . 



(A5) 



APPENDIX B: SUM FOR At > Ib 

When 6^1, one can simphfy in Eq. (|34|) 5* = dljdh, 
where 



oo oo 



e 

n=l n— 1 r— 



Using the integral representation of the F-function in 
Eq. (|Bip we have 



^-E 



r! r(-r/2) 



/ da;x-'"/2-i Ve-""- (B2) 



Note that the change of the order of sum and integral 
does not result in any extra terms.- The sum over n can 
now be trivially performed and using the relation 



1 



r(l + r/2) ,7rr, , , 

= ^ sin — , B3 

1 (— r/2) TT 2 



we get 



7 = _^/ (-1)'' ^ 

'ox _ 1) Z_. r(r + 1) I ./^ 



x(e--l) ^r(r + l) 



(B4) 



1=1 dx . 



(B5) 



This gives Xt in Eq. ([55]) . Alternatively, we could have 
expanded the logarithm term in and have kept only 
the first term in that expansion for 6^1; one has = 
—CTBI. Thus, one arrives at the same expression for 
XT as in Eq. 



For the saddle-point approximation of Fp [h) as defined 
in Eq. (|37p . we write: 



dxg{x)e-'''^'''^ 



g{xo)e' 



°' I dxexp 





h"{xo) 



{X - XqY 



(B6) 

with g{x) = [^/^{e^ — 1)] ^ and h{x) = b'^/4x -|-plogx, 
and 2:0 is defined by h'^xo) = (the prime refers to 
derivative). Simple manipulations following this scheme 
yield Eq. ([35]) and Eq. 
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